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1. Introduction 

Precision tests of the standard model from flavor factory and intensity frontier experiments 
can reveal new physics even if the accelerator energy isn't sufficient to create the new particles 
associated with the new physics. Here preliminary results of a lattice-QCD calculation of the 
nonzero-recoil form factor for the semileptonic process B — Dfv are presented. The principal goal 
is to determine \Vcb\ to a high precision. The theoretical uncertainty on \Vcb \ limits the precision of 
the unitarity triangle constraint from neutral kaon mixing. 

There are several methods for determining \Vcb\- (1) From the inclusive decay b — )• civ, pertur- 
bation theory and the operator product expansion provide an estimate of the inclusive decay rate, 
which is then combined with the measured decay rate to obtain \Vcb\ (2) From the exclusive 
semileptonic decay B — D*tv, lattice-QCD methods provide the hadronic contribution to the decay 
rate [Q]. For the special case of zero recoil, we have presented an unquenched calculation for this 
process three years ago with a preliminary update reported at the Lattice 2010 [@] and CKM 
2010 conferences (3) Likewise, from the exclusive semileptonic decay B — Div, lattice-QCD 
methods provide the hadronic contribution Here we report the first unquenched lattice-QCD 
calculation for this process at nonzero recoil. 

At present there is a 1.6 a disagreement between the value of \Vcb\ determined from inclusive 
decays and our recent preliminary value based on the decay B — D*£v at zero recoil, as illustrated 
in Fig. H Thus, further cross checks of and improvements in the theoretical calculations are needed. 
The value of \Vcb\ can be determined from the ratio of the measured decay rate to the calculated 
hadronic form factor at any chosen recoil energy. Typically, the uncertainty in the experimental 



measurement increases as the recoil energy vanishes [12], whereas the uncertainty in the lattice- 
QCD calculation decreases [|6[, ^ ^. Thus, to minimize the uncertainty in the ratio, it is best to 
combine results of the experimental measurement with those of the lattice-QCD calculation over 
the full range of available recoil energies. For this reason we undertake an analysis over a broad 
range. 

2. Methodology for determining \ Vcb\ 

The hadronic weak matrix element for this process is commonly parameterized as 

(D(/)|r^|B(/,)) = >i^[/i^-^°(w)(v + v')^+/j^-^^(w)(v-v');,] . (2.1) 
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Figure 1: The central values and errors resulting from three determinations of \Vcb\- (1) A "global fit" to 
the unitarity triangle including all inputs except \Vcb\ [^] [p^, (2) inclusive measurements summarized by 
the HFAG [ITTll, and (3) ourB ^ D* lattice-QCD calculation [F ~ 
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Here, w is velocity transfer, w = v • v'. The differential rate for the semileptonic decay B — )• Dfv is 



where ^b^d {w) is defined as 

<^b^d{w) = /^r°(w) - . (2.3) 

Thus, once we know and , we can determine the form factor ^^b^d(>v). Then, combin- 
ing it with the result of the decay rate from experiment, we will obtain \Vcb\- In order to construct 
and h^S^^ , we need the following quantities: 

h+{w) = /?+(p)[l -xy(p)/?_(p)], (2.4) 
h^{w) = /?+(p)[l -/?_(p)A/(p)], (2.5) 
/?+(p) = (D(p)|r4|B(0)), (2.6) 

_ (D(p)|r^|g(o)) 

^^P^ (D(p)|r4|B(0))' ^ ^ 

_ (D(p)ir^lD(O)) 
^^^P^ - (D(p)|r4|D(0)) ' ^^-^^ 

where y is the continuum hadronic weak vector current. 

In our calculation we use the local (nonconserved) lattice vector current {x) = {x)iy^^c {x) 
with i?{a) improved heavy quark fields *F/;(x) following [|^], and we renormalize it following a 
partly nonperturbative method, namely 



- Pv V Zybb^Vcc , (2.9) 



where the flavor-diagonal renormalization coefficients are computed nonperturbatively on the lat- 
tice via the conditions 

Zvbb{B\V^\B) = 1, (2.10) 
Zvcc{D\V^\D) = 1, (2.11) 

and Pv is computed perturbatively [0]. Because of cancellations among similar loop diagrams, we 
expect that pv is nearly equal to 1 |^. Preliminary results presented here omit the pv factor. 
At zero recoil (w = 1) we use the double ratio method |^: 

.,.|2_ ,2 {D\V'\B){B\V^\D) 

- {D\V^\D){B\V^\B) ■ ^^-^^^ 

This ratio suppresses a large part of the statistical fluctuations, and it builds in the current renor- 
malization ^ZybhZycc . 
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Table 1: Summary of all ensembles to be included in the full analysis. The ensembles in bold have been 
analyzed for this report. Light(strange) sea quark masses are denoted by m/(m/,). 

3. Results 

We report on results from an analysis of a large set of gauge field ensembles generated in the 



presence of 2 + 1 flavors of improved staggered (asqtad) quarks [ ]14[ , |15| ]. Further ensembles will 
be included in the future. Some key ensemble parameters are listed in Table |3[ 

On each ensemble we compute the three-point correlation functions relevant to the weak ma- 
trix elements and the two-point functions relevant to the propagation of the B and D mesons. For 
the heavy-light mesons we use heavy clover quarks in the Fermilab interpretation and light, im- 



proved staggered (asqtad) quarks [14, 15]. We use two types of interpolating operators for the 
B and D mesons, namely local and smeared using the IS Richardson wave function. We set the 
valence-quark masses equal to the sea-quark masses. The bare lattice charm and bottom quark 
masses are fixed by matching the kinetic masses of the D., and Bs mesons, respectively, to their 
experimental values. Some small adjustments will be needed to refine this tuning but are not yet 
included in this preliminary analysis. 

In terms of interpolating operators and ffi, for the B and D mesons and the vector current, 
the two-point and three-point functions are given by 

C2^^(f;p) =£exp(/p.x)(^(x,0^+(0)) 

= £4Z„(p)[exp(-£„(p)0+exp(-£„(p)(A^,-0)] (3.1) 

n 

4f^^°(?,r;p) = £exp(/p.y)(^z)(0)V^(y,0^1(x,r)) , (3.2) 

where 5„ = ±1 accounts for contributions that oscillate in t, and Nt is the lattice extent in the t 
dimension. Note that in the three-point function above, we have put the D meson at the origin and 
the B meson at (x, T). 



The matrix elements we want from Eq. ( [2.6| ) to Eq. ( [2.8| ) are obtained by factorization and 
reduction of the three-point functions. The reduction is based on the overlap coefficients Z„(p), 
energies and masses obtained in fits to the two-point functions. For example, the /?+(p) matrix 
element is obtained from 



V4 



I Zd{0)Ed ^Eot-\moT ^-mB{t-\T) (.3 3^ 

ZDip)mD 
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in the limit f — )• oo and T — f — > oo. In that limit it is a constant, /?_|_(p), independent of t. However, at 
finite t and T , complications from oscillating terms and excited states must be taken into account. 

Because they contain a light staggered quark, the meson interpolating operators excite even- 
as well as odd-parity channels. The unwanted even-parity states manifest themselves as terms that 
oscillate in t, corresponding to terms in the two-point and three-point functions with s„ = — 1. To 
suppress the effect of the oscillating terms in the three-point functions, we average in both t and T 
(after removing the dominant exponential dependence on t as in Eq. (^j|)) as follows: 



R=^-R{0,tJ)+^-R{0,t,T + \)+^-R{0,t + \J + l) . (3.4) 

Here, R refers to any of the quantities derived from three-point functions, such as /?+,/?_, and Xf. 

Having suppressed the contributions of the oscillating states to the three-point correlators, we 
must still account for contributions from excited states. For example, for D — 5- D, the leading 
contributions are 



x{D'\V4\D)^^y^Z;A^+ y^Z^^=={D\V4\D') VZ^, (3.5) 

V2wiz) v2£d \/2mD' 

where the primes indicate the overlaps and energies of the excited states of the same parity. We have 

neglected the doubly excited contribution, since the singly excited contribution is found already to 



be small. Thus, at finite t and T the ratio r;p) in Eq. (33) is not constant, but has a small 
contamination from excited states. We take this into account by fitting the resulting ratio to 

R+{p,t, T) = R+{p)exp{5mt) +Aexp(-A££,f) +5exp(-AmB {T - 1)) (3.6) 

as illustrated in Fig. ^ The fit parameter 5m should be zero, but it is introduced to allow for 
a small error in determining Ed and m^. The other terms involve AE^ = Eqi — Ed and Am^ = 
mg' —niB. These parameters are obtained from the fits of two-point functions. The results of two- 
point function fits then become prior central values and widths for fitting r;p). We use a 
similar method to obtain /?-(p) and xj(p). 

4. Chiral and continuum extrapolation 

To complete the analysis we will extrapolate in the light quark masses to their physical values 
and take the continuum limit, using staggered chiral perturbation theory [|T6|, ^] . At this prelim- 
inary stage, we don't show physical values and chiral fit in this report. The numerical data points 
for /j+ and h- are shown in Fig. |3[ We will need these data points to finish the extrapolation of 
physical value of and /j_ and the chiral fit. 

5. Future plans 

To complete the project, we need to (1) finish analyzing the available ensembles, (2) make 
adjustments for the tuning of the heavy quark masses, (3) compute the current renormalization (py 
factor), and (4) combine the experimental results for the differential decay rate with our results for 
the form factor to obtain | Vcb \ , using methods employed for | V,,^ | to extend w over the full kinematic 
range [p^. 
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Figure 2: A sample fit to r;p). Example three-point correlation fit for the 0.12 fm 0.14 m//m/, 
ensemble. The circles correspond to the local source, and diamonds, the smeared source. The fit result is 



obtained by a simultaneous fit to both correlators using Eq. ( |3.3| ). The horizontal lines represent the central 
value and error of the result for ^+(p). 
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Figure 3: Form factors /;+ (left) and h- (right) vs. w, omitting the matching factor pv, on the bold ensembles 
in Table |[ The color code is given in the legends. 
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